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Abstract. We use the adiabatic compression theory to build a physically well - motivated Milky Way mass model 
in agreement with the observational data. The visible mass of the Galaxy is distributed in a spheroidal bulge and 
a multi - components disc parametrized by three galactic parameters, the Sun distance to the galactic centre, Ro, 
the total bulge mass, Afbui ge , and the local disc surface density, E . To model the dark matter component, we 
adiabatically compress a Navarro, Frenk and White (NFW) halo (with concentration c and total mass M v i r ) for 
fixed values of the spin parameter, A, the fraction of the mass in baryons, m^, and the thin disc contribution to 
total angular momentum, ja- An iterative selection procedure is used to explore in very detail the wide space of 
parameters only selecting those combinations of {-Ro, Afbuige, E , A, mb, jd, c, M v i r } that give rise to a Milky Way 
model in agreement with the observational constraints. This analysis leads us to conclude that only models with 
Ro = 8.5 kpc, 0.8xl0 10 M < M bu i gc < 1.6xl0 10 M and 49 M pc" 2 < E < 56 M pc~ 2 can be reconciled 
with the set of observational constraints. As regard the parameters entering the adiabatic compression, we find 
0.03 < A < 0.10 and 0.04 < nib < 0.10, while the final estimates of the parameters describing the initial halo 
profile turn out to be 5 ~ c ~ 12 and 7xlO n M ~ M vir ~ 17x10" M Q (all at 95.7% CL). 
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1. Introduction 

The determination of the mass distribution of the Milky 
Way is a classical task of astronomy (jSchmidt 1 956 
ICaldwell fc Ostriker 198ll |Dehnen fc Binney 1998} . The 
usual framework for the origin of structures in 
the Universe is provided by the Cold Dark Matter 
(CDM) paradigm. This standard cosmological the- 
ory has proved very successful on large scale, in 
explaining both the abundance and clustering of 
galaxies UPeacock et al. 200ll IVerde et al. 2002fl and the 
power spectrum of the cosmic microwave background 
anisotropies l|de Bernardis et al. 2000(1 . but is experienc- 
ing a number of difficulties on the scales of galax- 
ies and dwarf galaxies. CDM paradigm predicts an 
over-abundance of satellites around the Milky Way 
and M31 by an order of magnitude ( |Klypin et al. 1999 
IMoore et al. 1999). While the presence of a photoion- 
izing background can solve this "sub-structure" prob- 
lem l|Somerville 200211 . difficulty about density profile 
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still remains. Dark matter dominated objects, such as 
dwarf and low surface brightness galaxies, show rotation 
curves inconsistent with the central density cusp pre- 
dicted by CDM cosmology ( |McCaugh fc de Block 19981 
lAlam et al. 2002|l . 

The understanding of galaxy formation also re- 
mains unsolved in the standard hierarchical model. 
Dissipationless CDM haloes are assumed to form bottom- 
up via gravitational amplification of initial density fluc- 
tuations. Gas carried with such haloes cools and con- 
tracts within them to form luminous, independent self- 
gravitating units, which can form stars, at the halo 
centres (|Fall fc Efstathiou 19801 IBlumenthal et al. 19861 
IMo et al. 1998(1 . The halo profile affects both gas cooling 
and infall since it determines the structural properties of 
the resultant heavy discs and dense nuclear bulges. Even 
if the growth of dark haloes is not much affected by the 
baryonic components, the halo gravitationally responds to 
the dissipative baryonic infall and the present day CDM 
haloes can have density profiles quite different from the 
original CDM prediction. 
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In this context, we want to address the density profile 
problem by fitting CDM models to the observed prop- 
erties of the Milky Way. The Milky Way seems to be a 
typical system on a mass scale of 10 12 M Q , mostly con- 
tributed by exotic particles, such as weakly interacting 
massive particles or axions. Direct searches for dark com- 
pact objects, such as MACHOs, in the Milky Way halo 
have been performed by the MACHO and EROS collab- 
orations through microlensing surveys. According to the 
MACHO group IjAlcock et al. 2000bjl . the most likely halo 
fraction in form of compact objects with a mass in the 
rang e 0.1-1 M is of ab out 20%; the EROS collabora- 
tion l|Lasserre et al. 2000|l has set a 95% confidence limit 
that objects less than 1 M Q contribute less than 40% of 
the dark halo. These upper limits on the fraction of com- 
pact objects are evidences favouring that Galaxy halo can 
be described by the standard cosmological approach. 

The very detailed data, obtained from many indepen- 
dent techniques, which characterize the Milky Way, make 
it a unique testing for theories of galactic structure and 
baryonic infall. Actually, this consideration has motivated 
Klypin et al. (2002) to perform an analysis similar to the 
one we present in this paper. Although conceptually anal- 
ogous to that of Klypin et al. (2002), the procedure we 
will use is much more detailed and allows to deeply ex- 
plore the whole parameter space. As a result, we will be 
able not only to investigate the viability of the standard 
CDM paradigm when applied to the Galaxy, but also to 
constrain a set of parameters that are, on the contrary, 
held fixed in the work of Klypin et al. (2002). 

The paper is organized as follows. In Sect. 2 we de- 
scribe the ingredients needed for our analysis and in Sect. 03 
we introduce briefly the adiabatic compression theory. 
Section 4 is devoted to the description of the observational 
constraints we have for the Milky Way, while in Sect. 5 we 
explain how we investigate the parameter space to select 
only models in agreement with the data. The impact of 
the different selection criteria on the parameter space is 
discussed in Sect. 6 and the results of our analysis are 
reported in Sect. 7. Section 8 is devoted to some final con- 
siderations. 

2. Model ingredients 

Traditional models of spiral galaxies usually include at 
least three dynamical components : a spheroidal bulge, an 
exponential disc and a dark halo. The lack of observational 
data precludes from using more detailed multi - component 
mass models since it is not possible to break the degener- 
acy among the many parameters involved. However, due 
to our privileged position, cinematic and photometric data 
on the Milky Way are numerous enough to require the use 
of a multi - component model in order to have a unified 
picture. To this aim, we model the mass distribution in 
the Milky Way introducing five components (the bulge, 
the thin and thick discs, the interstellar medium disc and 
the dark halo) that we describe in detail in the following 
subsections. 



2.1. The bulge 

The morphology of the Galactic bulge (defined as the 
spheroid within the galactic coordinates \l\ < 20° and |6| < 
10° region around the Galactic centre) is much harder to 
ascertain than that of the bulges in many external galax- 
ies, because of obscuration by interstellar dust due to our 
position in the Galactic plane. In the recent years, how- 
ever, striking images of the Galactic bulge have been ob- 
tained in the near infrared (at wavelengths of 1.25, 2.2, 
3.5 and 4.9 /im) by the DIRBE experiment on board the 
COBE satellite (jArendt et al. 19941 IWeiland et al. 1994jl 
allowing to study with good accuracy its structure. These 
images suggest that the stellar distribution in the bulge 
is bar -shaped, i.e. that the bulge is not rotationally 
symmetric. However, a bar -like structure can not be 
used since the standard adiabatic compression theory 
formally assumes that all the mass distributions have 
spherical symmetry. On the other hand, it is possible 
to achieve a relatively good agreement with the near in- 
frared photometric data also using spheroidal models for 
the bulge <|Kent et al. 19911 IDwek et al. 1995jl . Thus, we 
will describe the bulge as a spheroidal density distribution 
UDehnen fc Binney 19981 : 

where 

m 2 = (R 2 + z 2 /q 2 y/ 2 , (2) 

with R the galactocentric radius and z the height above 
the equatorial plane. Thus the density of the bulge is pro- 
portional to r~ 7 for r << ro, to r _/3 for r$ << r « r t 
and softly truncated at r = r t . Fitting of the model to 
the observed infrared photometric COBE/DIRBE data 
yields the values of four of the five bulge parameters 
UDehnen fc Binney 19981 : 

(3 = 7 = 1.8 , q = 0.6 , r = 1 kpc , r t = 1.9 kpc . 

The density normalization po is not determined from 
the fitting relation, but it is easily related to the total mass 
of the bulge, Mbuige- Integrating Eq.Q, one gets 

M h (m) = 0.518Af bu i g( V- 2 [0.6, 1.6, -0.27m 2 ] (3) 

where \F\ is the hypergeometric function and -Mbuige is 
related to the central density po (in M Q /pc 3 ) as : 

M bu i g e = 1.60851 xAwqpo . (4) 

In the adiabatic compression formalism, all the galaxy 
components are assumed to be spherical. To this aim, it 
is needed to describe the bulge mass distribution with a 
modified "spherical" version of Eq.®. A simple and rea- 
sonable way to solve this problem is to substitute Eq.J3J 
with the following one : 

A/rV) * M bulgc /(r) , (5) 
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f(r) = 1 - e~ l sb7r (l + 1.543r + 0.1898r 2 + 

0.6349r 3 -0.6109r 4 + 0.1491r 5 -0.01126r 6 ) , (6) 

where now r is the usual spherical radius. The spherical 
mass distribution in Eq.© has been denned so that the 
mean density inside the spherical radius r is the same as 
the one within the elliptical radius to, Mb (to), i.e. : 



K ph (r) 



M b (to) 



4/3 7t r 3 4/3 7r q to 3 

Using Eq. © instead of Eq. © introduces a systematic er- 
ror when describing the very inner regions of the Galaxy. 
Actually, this error is expected to be not a serious one 
since the bulge contribution to the total mass budget of 
the Galaxy is indeed small. Moreover, as we will see later, 
most of the observational constraints we will use probe a 
region of the Galaxy that is so far from the inner bulge 
dominated region that its dynamical effect could be de- 
scribed even by modelling this component as a pointlike 
mass. However, to further reduce this effect, wc will use 
Eq.© for the bulge mass distribution, but we consider 
the exact rotation curve ( |Binney fe Tremaine 1987| ) : 



r 1 

v 2 c (R) = iirGq / 
Jo 



Ph{m )m dm 
^/R 2 -(1-q^fr 



(J) 



To fully characterize the bulge we only need its to- 
tal mass, Mbuigc- Dwek et al. (1995) found Afbuigo = 
(1.3±0.5)xl0 10 M Q . 



2.2. The disc 

Contrary to the bulge, the structure of the Milky Way 
disc is quite easy to investigate given the large amount 
of available data. Disc is made up of three compo- 
nents, namely the thin and the thick stellar discs and 
the interstellar medium (ISM) disc. We model the two 
stellar discs with the usual double exponential profile 
HFreeman 19701 |Dehnen fe Binney 1998D : 



m \ Ed ( R l z l 
Pa{R,z) = — exp 

2z d \ R d z& 



(8) 



The total mass of each stellar disc is Aid = 27rEdi? 2 . For 
the adiabatic compression, we will adopt a 3-D mass dis- 
tribution such that : 



M d (r) = 2tt 



+ 00 



Pd{R, z)RdRdz 



(9) 



In Eq.©, Md(r) has spherical symmetry which is not for- 
mally correct since the disc is a highly flattened structure. 
Actually, this approximation introduces a negligible error 
as it is witnessed by the good agreement found between 
the predictions of the adiabatic compression theory and 
the numerical simulations IjJesseit et al. 2002)l . However, 
to evaluate the disc rotation curve and the vertical force 



we use the original flattened density distribution given in 
Eq.©. 1 

Although the simple formula in Eq.© may fit well 
the large-scale structure of the stellar discs, it is not able 
to reproduce the smaller scale density fluctuations, which 
are prevalent in the ISM component. Dame et al. (1987) 
have shown that there is very little interstellar matter be- 
tween the nuclear disc at 200 pc and the molecular ring 
at R = 4 — 5 kpc. These local fluctuations strongly affect 
the estimation of the Oort constants. Hence, we do not 
assume an analytical expression for the ISM disc density 
profile, but we use a third order polynomial interpolation 
of the data in Table Dl of Oiling & Merrifield (2001) also 
including a 23.8% helium contribution by mass. The 3-D 
mass distribution of the ISM disc is then evaluated as in 
Eq.©. The ISM disc rotation curve has been evaluated 
following the method described in Kochanek (2002). 

To fully characterize the disc model, we have to fix the 
geometrical parameters (Rd, Zd) and a value for the central 
(or the local) surface density. We fix the geometry of our 
stellar discs giving their scale - length and scale - height as 
follows ( |Dehnen fc Binney 1998| ) : 

thin disc : Rd = kRq kpc , z& = 180 pc ; 

thick disc : Rd = hRq kpc , Zd = 1000 pc . 

where Rq is the distance of the Sun to the galac- 
tic centre and K = (0.30±0.05) a scaling constant 
(Dc hnen fc Binney 1998| ). The exact value of i?o is quite 
uncertain, with most of the estimates ranging from 
7.0 to 8.5 kpc flKerr fc Lynden- Bell 1986| iReid 19931 
IQlling fc Merrifield 2000] ~ |Qlling fc Merrifield 200"T1 ). 
Given the importance of this parameter in the modeling 
of the Galaxy, we will explore models with different values 
ofi? - 

We use the following estimate for the local surface den- 
sity ( |Kuijken fe Gilmore 1989| ) : 

E = (48 ±8) Mq pc" 2 . 

Eq accounts for the contribution of the discs (both stellar 
and ISM), not taking into account the halo contribution. 
The total mass of each sub -disc is fixed specifying the 
fractional contribution of each one to Eq. Using the rela- 
tions given in Dehnen & Binney (1998), we fix: 



-Jthin 



14 E 



thick 



(Ro) 



while, for the ISM disc, the adopted value is 
jOlling fc Merrifield 20011 : 

SismORo) = 14.5 M Q pc" 2 . 

It is noteworthy that there are also other estimates 
of Eq significantly lower than the one used here 
l lUlling fc Merrifield 2001] IGerhard 207rl . 

1 To this aim we use a C++ code kindly provided us by 
W. Dehnen which implements a modified multipole technique 
developed by Kuijken & Dubinski (1994, see also Dehnen & 
Binney 1998). 
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2.3. The halo 

Observational data may be fitted by a wide range of mod- 
els, even unphysical ones. This is why there are a lot of 
different dark halo models which are claimed to describe 
well the density profile of this component. To obtain phys- 
ically interesting models, it is thus important to impose 
constraints based on a physical theory of halo formation 
and to select models which are both compatible with the 
data and also physically well motivated. From this point 
of view, numerical simulations of galaxy formation in hier- 
archical CDM scenarios are very helpful since they predict 
the initial shape of the dark matter distribution. In this 
paper, we assume a NFW profile (Nav arro et al. 1997}) as 
initial dark matter halo. The main properties of the NFW 
model are : 



p(r) 



x (1 + xf 



r/r s 



M(r) = 4 W 3 f(x) = M vi J(x)/f(c) , 
/(x)=ln(l+aO- 



1 + x ' 



r/r a 
AttS 



th 



^MPcrit^vir; 



(10) 

(11) 

(12) 
(13) 
(14) 



where c is the concentration parameter, M V i r the virial 
mass and r vir the virial radius 2 . The model is fully de- 
scribed by two independent parameters, which we assume 
to be c and M vir . A correlation between c and M vir has 
been found in numerical simulations ( Navar ro" et al. 19971 
IBullock et al. 20011 ICoiin et al. 2004)l . but we do not use 
such a relation since it is affected by a quite large scatter 
(-25%). 

The NFW model is not the only model proposed 
to fit the results of numerical simulations. Some au- 
thors IjMoore et al. 19981 |Ghigna et al. 2000| ) have pro- 
posed models with a central slope steeper than the NFW 
one. However, the difference between these models and 
the NFW one is very small for radii larger than 0.5%- 1% 
the virial radius and it is further washed out by the bary- 
onic infall. For these reasons, we will not consider models 
different from the NFW one. 



3. Adiabatic compression 

The present day dark matter halo has a different shape 
with respect to the original NFW model since the gravi- 
tational collapse of the baryonic matter, which forms both 
the bulge and the disc, changes the overall gravitational 
potential of the system. The halo structure is thus modi- 
fied by the forces of the collapsed baryons depending also 



2 The virial radius is defined such that the mean density 
within r v j r is <5th times the mean matter density of the universe 
p = S^MPcrit- We assume a flat universe with (Um,S1a, h) = 
(0.3,0.7,0.72) and S th = 337 flBryan fc Norman 1998^ . 



on the angular momentum of the galactic components 
(baryons and CDM). The effects of the baryonic infall 
are treated here following the approach of the adiabatic 
compression (Blumcnthal et al. 19861 IFlores et al. 1 993 
IMo et al. 1998(1 . If the disc is assembled slowly, we can 
assume that the halo responds adiabatically to the modifi- 
cations of the gravitational potential and it remains spher- 
ical while contracting. The angular momentum of the dark 
matter particles is then conserved and a particle which is 
initially at a mean radius r, ends up at a mean radius r 
where : 



M f (r) r = M(n) n 



(15) 



being M(ri) the initial total mass distribution within 
and Mf (r) the total final mass within r. Mf(r) is the sum 
of the dark matter inside the initial radius and the mass 
contributed by the baryonic components. We thus have : 



M f (r) = (l-m h )M(n) + . 



(16) 



where M{, ar (r) is the sum of the contributions from the 
baryonic components and is the fraction of the total 
mass that ends up in the visible components. Note that 
we are implicitly assuming that the baryons had initially 
the same mass distribution as the CDM particles and that 
those which do not form the luminous units still remain 
distributed as the CDM. 

For a given rotation curve, v c (R), the angular momen- 
tum of the thin disc is : 



J d = 2tt 



v c (R)Zf in (R)R 2 dR 



(17) 



The upper integration limit can be set to infinity since the 
disc surface density ~E^(R) drops exponentially and r vlr is 
much larger than the disc scale - length. We assume J<j to 
be a fraction jd of the initial angular momentum of the 
halo, J (i.e. Jd = jdJ)- J can be expressed in terms of a 
spin parameter A defined as : 



A = J|S| 1/2 G _1 M, 



-5/2 



with E, the total energy of the NFW halo, 



CM 2 

rp vir r 

~ Or ■ 
&l vir 

Some simple algebra allows us 
IMo et al. 199811: 



(18) 



(19) 



to rewrite Ea. I)17j) as 



Rd 



with : 



1 

71 



Ar vir / C 1/2 f R {X,c, m d , jd) 



l-l/(l + c) 2 -2(l + c)- 1 ln(l + c) 



[c/(l 



/i?(A,c, m d , jd) 



ln(l + c)] 2 



du 



(20) 



(21) 



(22) 
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In Eas. l|20l) and l|22l) . md is the fraction of the total mass 
competing to the thin disc, u = R/Rd and v V i r is the total 
circular velocity at the virial radius r vir . It is important to 
stress that the rotation curve entering Eq. I|17|) is the total 
one, i.e. it is : 



v 2 c (r) 



v 



c.bar 



(r) 



v. 



c,DM 



(r) , 



(23) 



where the first term is simply evaluated given the distri- 
bution of the baryonic components, while the latter is : 



vl IiM (r)=G[M i (r)-M hai (r)}/r . 



(24) 



The full set of equations allows one to determine the 
final distribution of the DM particles provided that a 
model for the density profiles of the baryons has been 
assigned and the parameters (A, mb, m<j, jd, c, M V i r ) have 
been fixed. We will see later in Sect. 4 that these param- 
eters are not all independent since it is possible to find 
some physically motivated relations among some of them. 

It is worth stressing that there is some debate about 
the validity of the adiabatic compression formalism. 
Jesseit et al. (2002) have found a substantial agreement 
between the final dark matter distribution in numerically 
simulated haloes and that predicted by the adiabatic com- 
pression approach. On the other hand, this result has been 
contradicted on the basis of a set of higher resolution nu- 
merical simulations recently carried out by Gnedin et al. 
(2004). According to these authors, the standard adia- 
batic compression formalism systematically overpredicts 
the dark matter density profile in the inner 5% of the virial 
radius. It is worth noting, however, that only one of the 
eight simulations considered by Gnedin et al. refers to a 
galactic (rather than a cluster) halo. The bottom panel of 
Fig. 4 in their paper shows that the adiabatic compression 
formalism overpredict the dark matter density less than 
~ 10% at r/r V j r ~ 0.1, while the error quickly decreases 
for larger values of r / r v ; r . This error is much smaller than 
the uncertainties we have on the observational quantities 
so that we are confident that using the standard adiabatic 
formalism does not introduce any bias in the results. 

4. Observational constraints 

We want to select mass models of the Milky Way that are 
physically well motivated and are in agreement with obser- 
vational data. We discuss in this section the observational 
constraints used to test each model. 

Only compact baryonic objects can cause microlens- 
ing events towards the Galactic bulge. The number 
of microlensing events observed towards the galactic 
bulge IjAlcock et al. 2000al |Popowski et al. 2000| ) deter- 
mines the minimum baryonic mass in the inner Galaxy 
that can yield the measured value of the optical depth 
t. For an axisymmetric mass density which decreases 
moving vertically away from the plane, it is possible 
to demonstrate that the minimum baryonic mass is 
UBinney fc Evans 200 ip : 



bar — 



G 



(25) 



where z — Ro tan b, with b the galactic latitude of the field 
target. This is only a lower estimate of the minimum bary- 
onic mass inside the solar circle and it is quite independent 
on the density profile of the baryonic components. Using 
the latest measured value of the optical depth towards the 
Baade window from the analysis of 52 events in which a 
clump giant is lensed, r = (2.0±0.4) x 10~ 6 (Popowski et 
al. 2000), we can thus select only those sets of parameters 
(Ro, -Mbuige, Eq) which produce a baryonic mass within 
i?o higher than what is predicted by Eq. 125(1 . Note that we 
are implicitly assuming that all the observed microlensing 
events are due to stellar lenses which is true only if there 
is no compact dark matter (such as MACHOs) in the disc. 
Should this assumption not to be true, the minimum bary- 
onic mass should be lowered. However, we stress that this 
constraint is not very selective so that changing the value 
of M 6 ™" does not affect our results. 

Whilst microlensing towards the bulge provides con- 
straints on the inner Galaxy, satellite dynamics and mod- 
eling of the Magellanic Clouds motion probe the Galaxy 
mass distribution on a large scale (~ 50- 100 kpc). Lin 
et al. (1995) have used the dynamics of the Magellanic 
Clouds to infer the mass of the Milky Way inside 100 kpc 
obtaining M(r < 100 kpc) = (5.5±1) x 10 11 M . This is in 
good agreement with the value found by Kochanek (1996) 
using escape velocity and motions of the satellite galax- 
ies which give (5 — 8)xl0 11 Mq. Taking into account the 
different techniques used and a dependence of the results 
on the halo modeling, we follow Dehnen & Binney (1998) 
assuming as our constraint : 



M(r < 100 kpc) = (7.0±2.5)xlO il M< 



(26) 



A third and quite efficient constraint is given by the 
Oort's constants defined as : 



.4 



v c (R) dv c (R) 
dR 



B = — 



R 



v c (R) dv c {R) 



R 



dR 



(27) 



(28) 



Dehnen & Binney (1998) have reviewed the estimates of 
the Oort's constants present in literature finally proposing 
the following values : 



A = (14.5±1.5) km s^kpc" 1 , 
B = (-12.5±2.0) km s^kpc" 1 , 
A — B = (27.0±1.5) km s^kpc" 1 



(29) 
(30) 
(31) 



From the definitions of the Oort's constants, it immedi- 
ately turns out that (^4 — B)Rq — v c (Rq). We can thus 
replace the constraints on A — B with a constraints on the 
local circular velocity : 



v c (R ) = i? o x(27.0±1.5) km s" 1 



(32) 



The vertical force K z at some height above the plane 
places a condition on the local mass distribution. Using 
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K stars as a tracer population, Kujiken & Gilmore (1989, 
1991) have deduced : 



K, 



\K Z (R ,1.1 kpc)| = 2ttGx(71±6) M pc~ 2 .(33) 



Formally, this estimate depends on the galactic constants 
R and v c (Rq), but Oiling & Merrifield (2001) have shown 
that the result is quite robust against variations in these 
parameters. Thus, we can use the previous value to con- 
strain our models. 

Finally, another constraint comes from the rotation 
curve of our Galaxy. We can reconstruct this quantity 
from the measurements of the velocity field. For an ax- 
isymmetric galaxy, the radial velocity relative to the local 
standard of rest, v r , of a circular orbiting object at galac- 
tic coordinates (/, b) and galactocentric radius R is related 
to the circular speed by : 



W(R) ee . Vr = ^v c (R) - v c (Ro) 
sin I cos o it 



(34) 



with the following relation between R and the distance d 
to the object : 



R 



d 2 cos 2 b - 



R 



2dRn cos b cos I 



Several studies are available with measurements of both d 
and v r for objects which ought to be on a nearly circular 
orbits, so that the Milky Way rotation curve v c (R) can be 
reconstructed. We use here the data on HII regions and 
molecular clouds in Brand & Blitz (1993) and the ones on 
a sample of classical Cepheids in the outer disc obtained 
by Pont et al. (1997). Following Dehnen & Binney (1998), 
we reject objects with either 155° < I < 205° or W < 
or d < 1 kpc when v r is very likely dominated by non - 
circular motions. 

5. Exploring the parameter space 

To apply the formalism of the adiabatic compression in or- 
der to have the present day halo mass profile, we need the 
initial halo shape and the today density profile of the bary- 
onic components. As discussed above, the baryons have 
been distributed in the bulge and in the three sub - discs 
according to the observational data. 

We have however still an indetermination on the bary- 
onic components since the Sun distance to the galactic 
centre Rq, the total mass of the bulge, Mb u i g e, the local 
surface density, £ , and the discs scale - lengths (fixed by 
the scaling constant k) are known with uncertainties. As 
a first step, we consider a grid of models with Rq ranging 
from 7.0 to 8.5 kpc in steps of 0.5 kpc, M bulge from 0.80 to 
1.80 xlO 10 Mq with a step of 0.125xl0 10 M , £ from 
40 to 56 Mq pc~ 2 in steps of 1.6 M Q pc~ 2 and k from 
0.25 to 0.35 in step of 0.05. Among these models we select 
only the ones which pass the test on the minimum bary- 
onic mass within Rq. We stress that changing the value 
of MV^ 1 to consider the (quite unlikely) possibility that 
some microlensing events are not due to stellar lenses have 



only a minor effect on this criterium. The lower is the value 
of M££ n , the higher is the number of models passing this 
preliminary test, but all of the models added by lowering 
M™ will be finally excluded by the selection procedure 
described later. 

The parameters (A, m\,, m<i, jd, c, M V i r ) characterize 
the dark halo and are needed to solve the adiabatic com- 
pression equations. Furthermore, they give scaling rela- 
tions between the baryonic components and the total mass 
distribution. First, we note that md can be expressed as 
function of m\, as : 



m d = 



— — — — xmb = Jd "»b , (35) 

M bu i ge + M thin + M thick + MiSM 



with Mthin, Mthick and Mism, respectively, the total mass 
of the thin, thick and ISM disc. It is also possible to ex- 
press the virial mass M v j r as function of M tn i n and m-b 
simply as l|Mo et al. 1998jl : 



M vir = 



M. 



thin 



M, 



thin 



m A f d m h 



(36) 



Equation (|36|) simply states that the final mass of the 
system is the same as the initial one. Thus, we are now re- 
duced to only four parameters, namely (A, m^, jd, c). We 
can further reduce the number of parameters observing 
that Eas. (|15(l . (|16|l and l|20(l are a system of three in- 
dependent equations which can be iteratively solved to 
determine the initial radius ri(r), the final mass Mf(r) 
and one of the four parameters, once fixed the remaining 
three ones. As we will see later, while there are some hints 
about the distribution of the other parameters, little is 
known about the value of the concentration c. We have 
thus decided to solve the set of equations with respect 
to (fj, Mf , c) having fixed the parameters (A, mb, jd)- The 
equations are highly non linear and must be solved itera- 
tively 3 , so that, to speed up the calculations, we have im- 
posed a priori that c should be in the range (5, 25) which 
is a quite conservative estimate for spiral galaxies similar 
to the Milky Way lj Jimenez et al. 2003|> . 

To explore in detail the space of parameters, we build, 
for each model with given values of (Ro, Mbuige, £©, k), 
a set of models individuated by the values of (A, mb, jd)- 
We briefly explain how we define the grid. We fix a value 
for the spin parameter A. The distribution of A for haloes 
generated in numerical N - body simulations is well ap- 
proximated by a log - normal distribution with param- 
eters nearly independent on the cosmological parame- 
ters, halo mass and redshift (Bar nes fc Efstathiou 1 987 
ILemson Kauffmann 19991 IVitvitska et al. 2002)) . Using 
the parameters in Vitvitska et al. (2001), the maximum 

3 Mo et al. (1998) proposed an approximate relation which 
can be used to estimate directly c from Eg. 1201 thus avoiding 
the iterative procedure. However, their approximation has been 
obtained neglecting the bulge and assuming a single exponen- 
tial disc instead of the three sub - discs we are using. We have 
checked that their formula may lead to underestimate strongly 
the value of c. 
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of the distribution is at A = 0.035 while there is a 90% 
probability that A is in the range (0.02,0.10). We thus 
let A change in this range in steps of 0.01. Next, we 
have to fix a value for m^. This parameter is poorly 
constrained since we only know that it cannot be larger 
than the universal baryon fraction ilb/^M- This latter has 
been inferred by observations involving completely differ- 
ent physical processes ((Turner 2002|l . The power spectrum 
of matter inhomogeneities from observations of large-scale 
structure is sensitive to fib/^M; the Two Degree Field 
Galaxy Redshift Survey has reported a value of 0.15±0.07 
IjPercival et al. 2001(1 . Measurements of the angular power 
spectrum of the CMBR also provide a very significant es- 
timate. The combined analysis in Jaffe et al. (2001) of 
several data set gives flb/^M = 0-186^q oos- We thus let 
mb change from 0.01 to 0.20 in steps of 0.01. Finally, 
we have to fix jd- This parameter also is not constrained 
either theoretically or by numerical simulations. On one 
hand, it is reasonable to assume jd = md and indeed this 
seems to be necessary to fit spiral galaxies rotation curves 
l|Mo et al. 1 998). On the other hand, numerical simula- 
tions have found jd/md significantly less than unity. We 
have varied this parameter from 0.005 to in steps of 
0.01. 

The grid we build in this way is quite detailed. For 
each given (Rq, Mb u i ge , k) the total number of sets 
(A, mb, jd) is ~ 800—1200 (depending on the value of nib), 
so that the parameter space is indeed checked intensively. 

To select among this large number of models we use, 
for each model parameterization, a multi-step procedure. 
First, we solve iteratively the set of Eas. (|15fl . H16fl and 
H20|) so that we have the full halo mass distribution. We 
are thus able to estimate Mf(r — 100 kpc), the total mass 
inside r = 100 kpc. A model passes to the next step only 
if this value is consistent with the estimate in Ea. (|26(l . 
In the second step, we compare the local circular velocity 
v c (Ro) to the constraint in Eq. l|32[l rejecting the model 
if there is no agreement. Next, the Oort's constant are 
evaluated and the model is retained if the values of A and 
B are in agreement with the values given by Eas. (|29|l and 
ijHUI) . Then, we compute K z> i.i and accept the model if the 
resulting value is in agreement with the estimate reported 
in Eq. (|33|l . For the surviving models, we estimate the x 2 
defined as : 



1 

N 



E 



u c,dat 



c, model 



(Ri) 



(37) 



where <r, is the error on the i - th measurements and the 
sum runs over the N data points. Note that we compute 
the x 2 as the last step of our selection procedure when all 
the model parameters are fixed so that we do not change 
their values in order to minimize the x 2 ■ We deem a model 
as acceptable if x 2 < 1-33. For our data set of N = 115 
entries, this corresponds to a 99.7% confidence level. 

Actually, as stressed by Oiling & Merrifield (2001), se- 
lecting among different models on the basis of the x 2 value 
is not a statistically correct procedure because the errors 



on both v c (R) and the Galactic constants, which enter the 
estimate of v c (R) through Ea. (|34l) . are not normal. This 
means that a high value of x 2 for a given model might be a 
consequence of an intrinsically wrong model or of the not 
normal origin of the errors. Furthermore, the data of Pont 
et al. (1997) on the radial velocities of the outer disc classi- 
cal Cepheids are given without any uncertainties, so that, 
for these data points, <7j is only determined by the propa- 
gation of the errors on the galactic constants and it is thus 
underestimated. So we have decided to still retain the x 2 ~ 
test a selection criterium, but we also consider the median 
statistics. As shown in Gott et al. (2001; see also Avelino 
et al. 2002, Chen & Ratra 2003, Sereno 2003), median 
statistics provide a powerful alternative to x 2 likelihood 
methods. Fewer assumptions about the data are needed. 
A proper median statistics assume that (i) experimental 
results are statistically independent; (ii) there are no sys- 
tematic effects. Statistical errors are not required to be 
either known or gaussianly distributed. Since our analysis 
is based on not normal errors, performing a test without 
using the errors themselves turns out to be a very con- 
servative approach. Furthermore, median statistics is also 
less vulnerable to the presence of bad data and outliers. 

To compute the likelihood of a particular set of pa- 
rameters, we count how many data points are above or 
below each model prediction and compute the binomial 
likelihoods. Given a binomial distribution, if we perform 
N measurements, the probability of obtaining k of them 
above the median is given by 



-N 



P(k) = 



Nl 



kl(N-k)\ 



(38) 



We count how many of the 115 experimental points are 
above the expected velocity rotation curve, and retain a 
model if the number of overestimates is between 43 and 72. 
Given the distribution in Eq. I|38[) . the probability that the 
median of 115 sorted entries falls in this range is 99.73%. 

It is worth to note that our final results do not depend 
on the sequence of the tests. 

6. Analysis of the selection criteria 

Before presenting the results of our analysis, it is inter- 
esting to investigate at what stage in our selection pro- 
cedure certain types of models are excluded, i.e. we want 
to study the impact of each criterium on the parameter 
space. To this aim, for a given set of galactic parameters 
(Ro, Afbuigo, Sq, k), we first select all the models with c in 
the range (5, 25) and then apply to this set of models the 
selection criteria introduced in the previous section sepa- 
rately. The main results of this analysis are presented in 
Figs. 1-5 and discussed below as regard the two param- 
eters entering the adiabatic compression equations, i.e. c 
and TOb- 

— Constraint on M{r < 100 kpc). As Fig. 1 shows, the 
application of this criterium tends to flatten the his- 
togram of the c values, being however slightly more 
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Fig. 1. Histogram of the number distribution of c and 
of rnb the for the models with (Rq, Mbuigc, ^0) n) = 
(8.5,1.0,54,0.3) and c between 5 and 25. Rq is in kpc, 
-Mbuige in 10 10 M Q , S© in M Q /pc 2 . The shaded histogram 
refers to the models passing the test on M(r < 100 kpc). 
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Fig. 2. Same as Fig. 1 for the test on v c (Ro). 



effective in cutting out models with values of c in 
the tails of the distribution. On the contrary, this cut 
tends to suppress the high end of the nib histogram, 
only retaining those models with smaller m^. This lat- 
ter result may be qualitatively explained considering 
Ea. H36|) which shows that the higher is mb, the lower 
is M v ; r and thus M(r < 100 kpc). 
Constraint on v c (Rq). Fig. 2 shows the impact of the 
test on the local circular velocity. High values of c 
are excluded by this constraint. As regard the mb his- 
togram, this selection criterium turns out to be almost 
orthogonal to the previous one since now low values of 
TOb are clearly disfavoured. 

Constraint on the Oort constants. This test turns out 
to be a sort of compromise between the two previous 
ones. On one hand, it is quite effective in excluding 
models with very low values of c (look at the lowest 
bin), as it is (with much less efficiency) for the con- 
straint on M(r < 100 kpc). On the other hand, it cuts 
away the low end of the mb histogram in the same 
way as the selection criterium based on the value of 
v c (Ro) do. This is a quite important result since it 
shows that the eventual exclusion of this constraint 
does not alter the final results of the multi-step selec- 
tion procedure. The presence of local fluctuations in 
the ISM disc density strongly affects the derivatives 
of the gravitational potential thus leading to possible 
errors in the evaluation of the Oort constants for a 
given set of galactic parameters. One could thus argue 
that it should be better to not use the Oort constants 
as a selection criterium, but only their difference and 
hence the local circular velocity. However, the use of a 
third order polynomial interpolation of the measured 
ISM disc density alleviates this problem so that we are 
confident that the estimated values of A, B are not cor- 
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Fig. 3. Same as Fig. 1 for the test on the Oort constants. 
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Fig. 4. Same as Fig. 1 for the test on K z .\.\. 

rupted. Furthermore, as we have observed comparing 
Fig. 3 with Figs. 1-2, applying only the constraint on 
A, B gives results that are consistent with those ob- 
tained using the two constraints on M(r < 100 kpc) 
and i> c (i?o) so that any systematic error in the estimate 
of the Oort constants is washed out in our multi-step 
procedure. 

— Constraint on K z ,i.i- This test works excluding models 
with values of c and in the tails of the distribution. 
This constraint is very effective when selecting among 
different galactic parameters leading to reject models 
with S Q < 49 M Q /pc 2 . Actually, the lower is S Q , the 
higher is the percentage of models excluded by the test 
on M(r < 100 kpc) or by that on the Oort constants. 
This is a clear evidence against models with low values 
of the local surface density in agreement with what the 
results of the application of the constraint on K Zt \_\ 
claim. 

— Constraint on y 2 . The application of this constraint 
allows to flatten the histogram of the c values lowering 
the peak in Fig. 5 for low values. However, the high 
end of the histogram is erased thus suggesting that 
models with very high values of c are not able to fit the 
Milky Way rotation curve. This is not an unexpected 
result since both N-body simulations and fitting of 
the adiabatically compressed NFW model to external 
galaxies show that small values of c are best suited to 
describe galactic dark haloes ( Jime nez et al. 2003)l . As 
concern , the constraint on the % 2 works as that on 
the local circular velocity selecting models with high 
values of this parameter. This is a reasonable result 
since both constraints are related to the same physical 
quantity, i.e. the rotation curve. The median statistics 
works the same way as the % 2 -test, but it is more 
stringent. 

As a final remark, we want to stress that the obser- 
vational data we have reviewed in Sect. 3 are "consensus" 
values with "consensus" errors so that these latter can not 
be treated as "statistical" uncertainties. This is the rea- 
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Fig. 5. Same as Fig. 1 for the test on x 2 value. 

son why we have decided to adopt a multi-step selection 
procedure instead of the usual x 2 minimization technique 
based on the definition of a x 2 entering all the constraints 
at the same time. Our filtering approach and the discus- 
sion presented in this section allow one to avoid all the 
problems connected with the statistics of not normal er- 
rors and makes it possible to understand how the results 
could vary changing one of the constraints. Actually, the 
procedure we have implemented is quite robust since the 
different constraints select different regions of the param- 
eter space thus allowing to narrow the ranges for both the 
adiabatic compression parameters and the galactic con- 
stants. 

7. Results 

The selection procedure we have employed is quite efficient 
allowing us to reject the most of the models. Since we have 
used two alternative test as last constraint (the x 2 value 
or the median statistics), we define two samples. Sample 
A contains those models passing all the selection criteria 
and having x 2 < 1-33, while Sample B is made out by the 
models passing the test on the median statistics. The final 
number of models is 116 in the Sample A and 34 in the 
Sample B from an initial set of ~ 10 6 . Actually, it turns 
out that the Sample B is a subset of Sample A, i.e. all 
the models in Sample B belong to Sample A too. This is 
expected since the median statistics is a more restrictive 
test than the x 2 analysis. 

The main results are shown in Figs. (6-9) and summa- 
rized in Tables 1 and 2. Note that the results obtained 
from the two samples are in perfect agreement so that we 
will not discuss them separately. We stress that the final 
samples of models may not be treated using the usual sta- 
tistical methods since our selection procedure is based on 
constraints on different observable that are "consensus" 
values with "consensus" errors. That is why we do not re- 
port as best estimate of the parameters their mean values, 
but the medians which is a more conservative approach. 
Because of this, the quoted 68% (95%) range must not be 
considered as the 1 — a (2 — a) confidence limit, but it is 
simply the range which the 68% (95%) of the values are 
within. With this caveat in mind, we discuss below the 
distribution of the various model parameters. 

1. Figure 6 shows the distribution of the values of the 
concentration parameter c. Small values are clearly 
favoured so that the median value for Sample A 
(Sample B) comes out to be c = 6.89 (6.44). It is worth 



Table 1. Median values, 68% and 95% regions for the 
model parameters of sample A. The total halo mass M v ; r 
is expressed in units of 10 11 Mq. 



Parameter 


Median 


68% range 


95% range 


c 


6.89 


5.49 - 9.45 


5.04 - 11.88 


M vir 


8.68 


7.53 - 10.20 


6.87 - 16.60 


m b 


0.08 


0.07 - 0.09 


0.04 - 0.10 


A 


0.06 


0.04 - 0.09 


0.03 - 0.10 


A' 


0.029 


0.026 - 0.035 


0.018 - 0.042 



Table 2. Median values, 68% and 95% regions for the 
model parameters of sample B. The total halo mass M v ; r 
is expressed in units of 10 11 Mq 



Parameter 


Median 


68% 


range 


95% 


range 


c 


6.48 


5.49 


- 7.26 


5.08 


- 9.75 


M v i r 


8.44 


7.58 


- 9.39 


7.03 


- 11.16 


m b 


0.08 


0.07 


- 0.09 


0.06 


- 0.09 


A 


0.06 


0.04 


- 0.09 


0.03 


- 0.10 


A' 


0.029 


0.027 


- 0.032 


0.023 


- 0.035 



17.5 
15 
12 . 5 

U 

10 

e 

3 7.5 
3 



2 . 5 






10 



12 



14 



Fig. 6. Histogram of the number distribution of the con- 
centration parameter c for the models in Sample A and 
Sample B (shaded histogram). 



stressing that the values of c we get are lower than ex- 
pected, but are not unrealistic. Jimenez et al. (2003) 
have fitted the rotation curves of 400 spiral galaxies by 
modelling them with an exponential disc and a dark 
halo obtained by adiabatically compressing the NFW 
profile. Their Fig. 1 shows the distribution of the con- 
centration c vs the total mass M v - a . There are indeed 
a lot of galaxies with values of c in the same range as 
the one found here. 

In Fig. 7, we plot the distribution of the values 
of the total mass M Y - a for the final set of mod- 
els. The median value for Sample A (Sample B) is 
8.68 (8.44) xlO 11 M , and the 95% confidence interval 
ranges from 5.04 (5.08) to 6.74 (9.75) xlO 11 Mq. These 
values are still in agreement with the most recent mea- 
surement of the Milky Way total mass estimated to 
be 1.9±^xl0 12 M Q (|Wilkinson & Evans 1999j> . It is 
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Fig. 8. Same as Fig. 6 for m^,. 



noteworthy, however, that the value given by Evans & 
Wilkinson is higher than usual being the most of the 
previous estimates lower than 10 12 Mq (see, e.g., Fig. 3 
in Zaritsky 1999). Sample A also contains few models 
with very high values of M v - n -. The value of M v - lr we 
have obtained may thus be considered quite reason- 
able. It is also possible to see that a clear correlation 
exists between c and M V j r , with the lower values of the 
concentration parameter c giving rise to the highest 
values of the total mass. This correlation may reflect 
the existence of some covariance among the model pa- 
rameters in the analysis procedure. 
The distribution of is shown in Fig. 8. The 
median value is 0.08 for both Sample A and 
B, while the 95% range turns out to be 0.07- 
0.09. According to some authors (JMo et al. 1998 
I Jimenez et al. 20031 |Klypin et al. 20 02), should be 
written as exOb/^M with e indicating the efficiency 
of the transfer of baryons from the initial halo to the 
bulge and disc. If this were correct, our results should 
mean e ~ 37 - 48% for n h /n M = 0.186. The efficiency 
is not one and this rises the problem of understanding 
where the missing baryons are. 

The median value of the spin parameter A is 0.06 for 
both Sample A and B and the 95% range is 0.04 — 

0. 09. Models with low spin are erased by the selection 
procedure. We do not report any result for j<j since this 
parameter is found to be degenerate with A and rrid, 

1. e., for fixed values of (mt,, -Ro, Mbuige, X©, k), models 
having the same value of A' = Xjd/md give rise to 
the same present day halo mass profile. This result is 
not unexpected since in the only equation containing A 
and jd, Eq. (|20|) . these two parameters appear only in 
A'. The median value for this latter is A' = 0.029 with 
a 95% range going from 0.018 to 0.042 for Sample A 
and from 0.023 to 0.035 for Sample B. 

Figure 9 shows the histogram of the number of mod- 
els according to the values of Mb u i ge and E Q . We can 
draw some interesting limits. First, we observe that 
we are not able to significantly constraint the bulge 
total mass. Indeed, the median value of the distribu- 
tion turns out to be Mb u i°-c = l.OxlO 10 Mq, but the 
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Fig. 9. Histograms of the number distribution of models 
according to the values of Mb u i ge (10 10 Mq, left panel) 
and Sq (Mq pc~ 2 , right panel). Shaded regions refer to 
Sample B. 

95% range for Sample A does not reject any model. It 
is worth noting, however, that values of Mbuige in the 
range 0.8— 1.4xl0 10 Mq seems to be favoured. This is 
confirmed by analysing Sample B that excludes models 
with Mbuigo > 1.4xl0 10 Mq. A similar analysis leads 
to the following constraints on the disc local surface 
density : 

49 Mq pc -2 < E Q < 56 Mq pc~ 2 

with 54 Mq pc -2 as median value. If we subtract the 
contribution of the ISM disc surface density, we get for 
the stellar disc X* ~ 40 Mq pc~ 2 in agreement with 
the consensus value (X* = 35±10 Mq pc~ 2 ) proposed 
by Oiling & Merrifield (2001). Only 6 out of 116 models 
in Sample A have i?o = 8.0 kpc, while in all the other 
cases it is Rq = 8.5 kpc so that we may safely consider 
i?o = 8.5 kpc as our final estimate of this galactic con- 
stant. This conclusion is strengthened by noting that 
all the models in Sample B have Rq = 8.5 kpc. This 
value is somewhat higher than expected since most re- 
cent estimates predict values in the range 7.0 - 8.0 
kpc IjReid 19931 |Qlling fc Merrifield 2000} . However, 
the estimate of Ro is somewhat model dependent since 
its value is linked to the halo flattening. For instance, 
Oiling & Merrifield (2001) have found that it is pos- 
sible to build galaxy models with Ro ~ 7.0 kpc if 
the halo is close to spherical which indeed is our 
case. Finally, there are no models with n ^ 0.30 
so that we conclude that the discs scale -length is 
Rd = 0.30xi? = 2.55 kpc. This is lower than the fidu- 
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Fig. 10. Rotation velocity curve in the Galactic 
plane for the models {£©, Mb u i ge) c, mb, jd, A} = 
{54.4M pc~ 2 , 1.05x 10 10 M Q , 5.55, 0.08, 0.035, 0.04} 
(full line) and {E Q , Af b uigc, c, m b , jd, A} = 
{56.OM pc- 2 , 9.25xl0 9 M Q , 6.56, 0.09, 0.045, 0.04} 
(dashed line). 

cial value (3.5 kpc) suggested by Binney & Tremaine 
(1987) and often adopted in many Milky Way disc 
modelling. However, we note that a value of ~ 2.5 is 
favoured by both star count models l|Robin et al. 1992 
Orti z fc Lepine 1993| ) and integrated NIR luminosity 
profile IjFreudenreich 19981 |Binney et al. 1997| ). 

Finally, we wish to comment again on the use of the 
X 2 as a selection criterium. Actually, it is difficult to say 
whether a certain value of \ 2 for a given model is due to 
some systematic error (as having neglected the triaxial 
structure of the bulge and the flattening of the halo) or to 
a problem with the parameter values or to the not normal 
nature of the estimated errors on the rotation curve. 
However, the agreement among the results obtained 
using both the \ 2 analysis and the median statistics 
(that turns out to be a more selective criterium) is a 
convincing evidence that the excluded models have been 
rejected as a consequence of physical problems, either due 
to the need for a more careful description of the inner 
Galaxy or to an intrinsically wrong combinations of the 
parameters (A,mb,jd,c). This makes us confident in the 
results. As an example, in Fig. ^J] we plot the rotation 
curve in the Galactic plane for two models, the one with 
the lower value of the % 2 ({£©, -Mbuigc, c, mb, jd, A} = 
{54.4M oPC - 2 , 1.05x 10 10 M Q , 5.55, 0.08, 0.035, 0.04}) 
and the one with median test value 
equal to 58 ({£©, M bulg0 , c, m b , j d , A} = 
{56.OM pc- 2 , 9.25xlO 9 M , 6.56, 0.09, 0.045, 0.04}). 
As we can see, the difference in the rotation curve 
between the two models is really negligible. We remind 
that each set of model with passes the selection procedure 
has the same statistical weight. 

8. Conclusions 

Modelling Milky Way is one of the classical tasks of as- 
tronomy. In this paper, we have applied the adiabatic com- 



pression method in the framework of the CDM structure 
formation to build models of the Galaxy in agreement 
with the observational data and, at the same time, well 
motivated by numerical simulations of galaxy formation. 
Exploring in very detail the parameter space, we have fi- 
nally selected a set of models which has allowed us to 
draw some interesting constraints on (i) the galactic pa- 
rameters, namely the Sun distance to the Galactic centre 
Ro, the total mass of the bulge Mb u i go , the disc local sur- 
face density £© and the disc scale -length Rd, (ii) the pa- 
rameters entering the adiabatic compression, i.e. the spin 
parameter A, the fraction of the mass in baryons m b and 
the thin disc contribution jd to the total angular momen- 
tum and (in) the parameters describing the initial halo 
profile (i.e. c and M V i r )- 

It should be interesting to compare our results with 
previous ones in literature, but a straight comparison is 
not possible because of the different approaches followed 
in the analyses. Klypin et al. (2002) have used the adi- 
abatic compression theory to build Milky Way models 
which are in agreement with observational data starting 
from a NFW halo. Even if the final aim is the same as our 
one, the approach used is radically different since they 
give a value for the concentration of the initial NFW halo 
and then find by trial and errors the values of the disc 
and bulge parameters so that the model is in agreement 
with the data. In our approach, the visible components 
are fixed from the beginning and the halo parameters 
are determined later. Klypin et al. (2002) finally examine 
only four models. Their favoured one is radically differ- 
ent from our ones since it has a longer disc scale - length 
(3.5 vs 2.55 kpc) and a higher concentration (12 vs. 6.48). 
However, it is noteworthy that their Model A2 seems to 
be more in agreement with our results, having c = 5 and 
M v i r = 7.1 xlO 11 M Q , and the authors find that also this 
model is able to fit the data. We stress, however, that 
most of the disagreements between our results and those 
of Klypin et al. (2002) may derive from a different way 
of selecting the models compatible with the observational 
constraints. 

A somewhat surprising result of our analysis is the 
quite small values of the concentration parameter c that 
turns out to be smaller than 11.88 (9.75) at 95% for 
Sample A (Sample B). Using the relation between M v ; r 
and c in Bullock et al. (2001), for values of M v - lr in the 
95% range reported in Table 1, one should expect c in 
the range 14.78 - 16.04. Considering a ~ 25% scatter 
( Col in et al. 20041 , our median c is more than 2a smaller 
than the above lower limit. Moreover, Fig. 6 shows that 
none of the models have such high values of c so that we 
may conclude that the quoted relation c — M v i T found in 
numerical simulations is not verified by our observation- 
ally selected Milky Way mass models. We have verified 
that this result is independent on our choice of cosmolog- 
ical parameters by repeating analysis for an Einstein-de 
Sitter model of universe (J7m = E = 0). 

Although surprising, this result is not fully unex- 
pected. Fitting adiabatically compressed NFW models 
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to the rotation curve of a sample of 400 spiral galaxies, 
Jimenez et al. (2003) have determined the values of c and 
M v ; r for these haloes. Looking at their Fig. 1, one sees that 
there are a lot of galaxies with values of M V - 1T in the range 
determined by us for the Milky Way. For these galaxies, 
the concentration parameter turns out to be of the same 
order as those obtained here so that the relation c — A/ V j r 
is not satisfied for these galaxies too. We are thus confi- 
dent that the disagreement we have found for the selected 
Milky Way models is not a result of the our procedure, but 
a possible shortcoming either of the numerical simulations 
or of the ACDM paradigm itself. 

The procedure presented may be extended to other 
spiral galaxies. A possible target is the Andromeda galaxy 
whose visible components can be modelled in detail thanks 
to the available photometric and cinematic data. It is 
also interesting to apply our procedure to low surface 
brightness (LSB) galaxies. Modelling in detail the visi- 
ble components in these galaxies is much difficult because 
of the paucity of the data, but these systems are prob- 
ably dark matter dominated so that uncertainties in the 
baryonic components should not affect systematically the 
main results. Moreover, combining the constraints from 
many galaxies will allow reducing the systematics con- 
nected with the disc modelling. To this aim, a sample of 
LSB with high resolution rotation curve should be ideal 
since it allows to decouple the disc contribution from the 
halo one l|de Blok fc Bosma 2002) . 

The analysis carried out may be further refined in 
order to investigate whether the constraints on the pa- 
rameters are affected by some simplifying hypotheses. We 
have implicitly assumed that there is no exchange of an- 
gular momentum between the dark matter particles and 
the baryons infalling in the disc and bulge. An approxi- 
mate analytical approach to this problem has been devel- 
oped ( |Klypin et al. 2002] ) and it should be thus interesting 
to repeat our analysis by including this effect. However, 
the details of how the baryons are transferred from the 
halo to the disc depend on the physics of star formation 
and supernovae explosion and some kind of energy feed- 
back is needed to eliminate the so - called cooling catas- 
trophe jBalogh et al. 2001] ). In order to investigate these 
effects, numerical simulations are needed, but the physics 
of the process may be also described with semi - analytical 
modelling IjCole et al. 20001 Ivan den Bosch 2002jl . These 
methods are best suited to be included in our procedure, 
thus making it possible to put stronger constraints both 
on the adiabatic compression parameters (A,mb,Jd) and 
on the initial NFW halo ones (c, M v - K ). 

We would like to conclude with a general considera- 
tion. The aim of this paper has been to build a Milky 
Way mass model in agreement with the observational data 
and motivated by some physical background. The simple 
adiabatic compression formalism has allowed us to take 
into account analytically the effect of baryonic infall in 
a way which is consistent with more complex numerical 
simulations. We have finally found that such a model in- 
deed exists although there is a disagreement with the pre- 



dicted concentration parameter. In our opinion, this could 
be considered as an indirect evidence suggesting that the 
hierarchical CDM scenario of galaxy formation is essen- 
tially correct since its predictions are in agreement with 
the data on the best studied galaxy, the Milky Way. We 
are thus confident that solutions to the problems of this 
cosmological model do not require modifications of dark 
matter properties, but have to be searched in astrophysi- 
cal phenomena. 

Acknowledgements. We acknowlege an anonymous referee for 
having pointed out a major error in our procedure and for 
his comments that have helped us to significantly improve the 
paper. We warmly thank Walter Dehnen for having furnished 
us his code for computing the gravitational potential of the 
thin and thick discs. It is a pleasure to thank also Shude Mao 
and Yogesh Joshi for some useful comments. 

References 

Alam, S.M.K., Bullock, J.S., Weinberg, D.H., 2002, ApJ, 572, 
34 

Alcock, Ch., Allsman, R.A., Alves, D.R., Axelrod, T.S., Becker, 
A.C., Bennett, D.P., Cook, K.H., Drake, A.J. et al.2000a, 
ApJ, 541, 734 

Alcock, Ch., Allsman, R.A., Alves, D.R., Axelrod, T.S., Becker, 
A.C., Bennett, D.P., Cook, K.H., Dalai, N. et al.2000b, 
ApJ, 542, 281 

Arendt, R.G., Berriman, G.B., Boggess, N., Dwek, E., Hauser, 
M.G., Kelsall, T., Moseley, S.H., Murdock, T.L. et al. 1994, 
ApJ, 425, L85 

Avelino, P.R, Martins, C.J.A.R, Pinto, P., 2002, ApJ, 575, 989 
Balogh, M.L., Pearce, F.R., Bower, R.G., Kay, S.T. 2001, 

MNRAS, 326, 1228 
Barnes, J., Efstathiou, G. 1987, ApJ, 519, 375 
Binney, J., Evans, N.W. 2001, MNRAS, 327, L27 
Binney, J., Gerhard, O.E, Spergel, D.N. 1997, MNRAS, 288, 

365 

Binney, J., Tremaine, S. 1987, Galactic dynamics, Princeton 

University Press 
Blumenthal, G.R., Faber, S.M., Flores, R., Primack, J.R., ApJ, 

301, 27 

Brand, J., Blitz, L. 1993, A&A, 275, 67 

Bullock, J.S., Kolatt, T.S., Sigad, Y., Somerville, R.S., 
Kravtsov, A.V., Klypin, A., Primack, J.P., Dekel, A. 2001, 
MNRAS, 321, 559 

Bryan, G.L., Norman, M.L., 1998, ApJ 495, 80 

Caldwell, J.A.R., Ostriker, J. P., 1981, ApJ, 251, 61 

Chen, C, Ratra, B., 2003, PASP, 115, 1143 

Cole, S., Lacey, C.G., Baugh, CM., Frenk, C.S. 2000, MNRAS, 
319, 168 

Colin, P., Klypin, A., Valenzuela, O., Gottlober, S. 2004, ApJ, 
612, 50 

de Bernardis, P., Ade, P.A.R., Bock, J. J., Bond, J. R., Borrill, 
J., Boscaleri, A., Coble, K., Grill, B.P., 2000, et al., Nature, 
404, 955 

Dame, T.M., Ungerechts, H., Cohen, R.S., de Geus, E.J., 
Grenier, I.A., May, J., Murphy, D.C., Nyman, L.A., 
Thaddeus, P. 1987, ApJ, 322, 706 
de Block, W.J.G., Bosma, A. 2002, A&A, 385, 816 
Dehnen, W., Binney, J. 1998, MNRAS, 294, 429 



V. F. Cardone and M. Sereno: Modeling the Milky Way 



13 



Dwek, E., Arendt, R.G., Hauser, M.G., Kelsall, T., Lisse, CM., 
Moseley, S.H., Silverberg, R.F., Sodroski, T.J., Weiland, 
J.L. 1995, ApJ, 445, 716 

Fall, S.M., Efstathiou, G., 1980, MNRAS, 193, 189 

Flores, R., Primack, J.R., Blumenthal, G.R., Faber, S.M. 1993, 
ApJ, 412, 443 

Freeman, K.C 1970, ApJ, 160, 811 

Freudenreich, H.T. 1998, ApJ, 492, 495 

Gerhard, O. 2002, Mass distribution in our Galaxy, in Matter- 
in the Universe, Jetzer Ph. et al. eds. 

Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., 
Stadel, J. 2000, ApJ, 544, 616 

Gnedin, O.Y., Kravtsov, A.V., Klypin A.A., Nagai, D. 2004, 
ApJ, 616, 16 

Gott III, J.R., Vogeley, M.S., Podariu, S., Ratra, B., 2001, ApJ, 
549, 1 

Jaffe, A.H., Ade, P.A.R., Balbi, A., Bock, J.J., Bond, J.R., 
Borrill, J., Boscaleri, A., Coble, K., et al., 2001, Phys. Rev. 
Lett., 86, 3475 
Jesseit, R., Naab, T., Burkert, A. 2002, ApJ, 571, L89 
Jimenez, R., Verde, L., Oh, S.P. 2003, MNRAS, 339, 243 
Kent, S.M., Dame, T.M., Fazio, G. 1991, ApJ, 378, 131 
Kerr, F.J., Lynden-Bell, D. 1986, MNRAS, 221, 1023 
Klypin, A., Kravtsov, A.V., Bullock, J.S., Primack, J.R. 2001, 
ApJ, 554, 903 

Klypin, A., Kravtsov, A.V., Valenzuela, O., Prada, F., 1999, 
ApJ, 522, 82 

Klypin, A., Zhao, H.S., Somerville, R.S. 2002, ApJ, 573, 597 
Kochanek, C.S. 1996, ApJ, 457, 228 
Kochanek, C.S. 2002, astro - ph/0108162 
Kuijken, K., Dubinski, J. 1994, MNRAS, 269, 13 
Kuijken, K., Gilmore, G. 1989, MNRAS, 239, 605 
Kuijken, K., Gilmore, G. 1991, ApJ, 367, L9 
Lasserre, T., et al., 2000, A&A, 355, L39 
Lemson, G., Kauffmann, G. 1999, MNRAS, 302, 111 
Lin, D.N.C., Jones, B.F., Klemola, A.R. 1995, ApJ, 439, 652 
McGaugh, S., de Block, W.J.B., 1998, ApJ, 499, 41 
Mo, H.J., Mao, S., White, S.D.M., 1998, MNRAS, 295, 319 
Moore, B., Governato, F., Quinn, T., Stadel, J. 1998, ApJ, 499, 
L5 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., 

Stadel, J., Tozzi, P., 1999, ApJ, 524, L19 
Navarro, J.F., Frenk, C.S., White, S.D.M. 1997, ApJ, 490, 493 
Oiling, R.P., Merrifield, M.R. 2000, MNRAS, 311, 361 
Oiling, R.P., Merrifield, M.R. 2001, MNRAS, 326, 164 
Ortiz, R., Lepine, J.R.D. 1993, A&A, 279, 90 
Peacock, J.A., et al., 2001, Nature, 410, 169 
Percival, W.J., Baugh, CM., Bland-Hawthorn, J., Bridges, T., 

Cannon, R., Cole, S., Colless, M., Collins, C, et al., 2001, 

MNRAS, 327, 1297 
Pont, F., Queloz, D., Bratschi, P., Mayor, M. 1997, A&A, 318, 

416 

Popowski, P. et al. 2000, astro - ph/0005466, in Microlensmg 
2000: A New Era of Microlensing Astrophysics, ASP Conf. 
Proa, Vol. 239, eds. J.W. Menzies and P.D. Sackett, San 
Francisco (USA) 

Reid, M.J. 1993, ARA&A, 31, 345 

Robin, A.C, Creze, M., Mohan, V. 1992, A&A, 265, 32 

Sereno, M., 2003, A&A, 412, 341 

Schmidt, M., 1956, Bull. Astr. Inst. Neth., 13, 15 

Sommerville, R., 2002, ApJ, 572, L23 

Turner, M.S., 2002, ApJ, 576, L101 

van den Bosch, F.C. 2002, MNRAS, 332, 456 



Verde, L., Heavens, A.F., Percival, W.J., Matarrese, S., Baugh, 

CM., Blass - Hawthorn J., Bridges, T., Cannon, R. et al., 

2002, MNRAS, 335, 432 
Vitvitska, M., Klypin, A., Kravtsov, A.V., Wechsler, R.H., 

Primack, J.R., Bullock, J.S. 2002, ApJ, 581, 799 
Weiland, J.L., Arendt, R.G., Berriman, G.B., Dwek, E., 

Freundenreich, H.T., Hauser, M.G., Kelsall, T., Lisse, CM. 

et al. 1994, ApJ, 425, L81 
Wilkinson, M.I., Evans, NW. 1999, MNRAS, 310, 645 
Zaritsky, D. 1999, in the Galactic halo, eds. B.K. Gibbon, T.S. 

Axelrod and M.E. Putnam, San Francisco 



